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Abstract 

Let k be a field. Then Gaussian elimination over k and the Eu- 
clidean division algorithm for the univariate polynomial ring k[x] allow 
us to write any matrix in SL n (k) or SL n (k[x]), n > 2, as a product of 
elementary matrices. Suslin's stability theorem states that the same 
is true for the multivariate polynomial ring SL n (k[xi, ■ ■ ■ ,x m ]) with 
n > 3. As Gaussian elimination gives us an algorithmic way of finding 
an explicit factorization of the given matrix into elementary matrices 
over a field, we develop a similar algorithm over polynomial rings. 



1 Introduction 

Immediately after proving the famous Serre 's Conjecture (the Quillen-Suslin 
theorem, nowadays) in 1976 [|TTJ, A. Suslin went on |12|| to prove the following 



.fTi-analogue of Serre 's Conjecture which is now known as Suslin's stability 
theorem: 

Let R be a commutative Noetherian ring and n > max(3, dim(R) + 
2). Then, any n x n matrix A = (/y) of determinant 1, with /y 
being elements of the polynomial ring R[x\, . . . , x m ], can be writ- 
ten product of elementary matrices over R[ 
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Definition 1 For any ring R, an n x n elementary matrix Eij(a) over R is 
a matrix of the form I + a - where i ^ j,a G R and is the n x n matrix 
whose component is 1 and all other components are zero. 

For a ring R, let SL n (R) be be the group of all the n x n matrices of deter- 
minant 1 whose entries are elements of R, and let E n (R) be the subgroup 
of SL n (R) generated by the elementary matrices. Then Suslin's stability 
theorem can be expressed as 

SL n (R[xi, . . . , x m ]) = E n (R[xi, . . . , x m ]) for all n > max(3, dim(i?) + 2). (1) 

In this paper, we develop an algorithmic proof of the above assertion over a 
field k. By implementing this algorithm, for a given A e SL n (k[xi, . . . , x m ]) 
with n > 3, we are able to find those elementary matrices Ei,...,E t e 
E n (k[xi, . . . , x m ]) such that A — E\- • ■ E t . 

Remark 1 If a matrix A can be written as a product of elementary matrices, 
we will say A is realizable. 



• In section 2, an algorithmic proof of the normality of E n (k[xi, . . . , x m ]) 
in SL n (k[xi, . . . , x m ]) for n > 3 is given, which will be used in the rest 
of paper. 

• In section 3, we develop an algorithm for the Quillen Induction Process, 
a standard way of reducing a given problem over a ring to an easier 
problem over a local ring. Using this Quillen Induction Algorithm, we 
reduce our realization problem over the polynomial ring R[X] to one 
over Rm[X\s, where R = k[x\, . . . , x m -i] and M is a maximal ideal of 
R. 

• In section 4, an algorithmic proof of the Elementary Column Property, 
a stronger version of the Unimodular Column Property, is given, and 
we note that this algorithm gives another constructive proof of the 
Quillen- Suslin theorem. Using the Elementary Column Property, we 
show that a realization algorithm for SL n (k[x 1 , . . . , x m ]) is obtained 
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from a realization algorithm for the matrices of the following special 
form: 

/ p q \ 

r s G SL 3 (k[x 1 , . . .,x m ]), 
\0 lj 

where p is monic in the last variable x m . 

In section 5, in view of the results in the preceding two sections, we note 
that a realization algorithm over k[xi, . . . be obtained from a 

/ p q 

realization algorithm for the matrices of the special form r s 

\0 1 

over where R is now a local ring and p is monic in X. A realiza- 

tion algorithm for this case was already found by M.P. Murthy in [f|]. 
We reproduce Murthy 's Algorithm in this section. 

In section 6, we suggest using the Steinberg relations from algebraic 
.fT-theory to lower the number of elementary matrix factors in a factor- 
ization produced by our algorithm. We also mention an ongoing effort 
of using our algorithm in Signal Processing. 



2 Normality of E n (k[x\, . . . , x m ]) in SL n (k[xi, . . . , x m ] 



1 + xy x 
—y 2 1 — xy 



is not realizable, but 



Lemma 1 The Cohn matrix A = 
^ J^J G SL 3 (k[x,y]) is. 

Proof: The nonrealizability of A is proved in and a complete algorithmic 
criterion for the realizability of matrices in SL 2 (k[xi, . . . ,x m ]) is developed 
13 ]. Now consider 



A 
1 



/ 1 + xy x 2 
—y 2 1 — xy 

V o oi 



(2) 



3 



1 + xy x 2 \ fx 
Noting that | — y 2 1 — xy = / + — y | • (y, x, 0), we see that the 

1/ V 

realizability of this matrix is a special case of the following Lemma 3. □ 



Definition 2 Let n > 2. A Cohn-type matrix is a matrix of the form 

I + av • (fj-ei — Wjej) 

E (k[x!, . . . ,x m ]) n , i < j G {!,..., n}, a e fc[xi, . . .,x m ], 



where v 



and e ; = (0, . . . , 0, 1, 0, . . . , 0) with 1 occurring only at the i-th position. 

Lemma 2 Any Cohn-type matrix for n > 3 is realizable. 
Proof: First, let's consider the case % — \,j — 2. In this case, 

■ (^,-^1,0, ...,0) 
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J + a 

/ 1 + at> !f 2 



-av\ 



0\ 



1 — aviv 2 • • • 



af 3 W2 -av 3 vi 

\ av n v 2 -av n vi 

f 1 + at> if 2 — a^i 

av 2 1 — av iv 2 



ra-2 



/ 

o\ 













'n-2 



V 



o o 

So, it's enough to show that 



Y[ Ei^av^E^-aviVx)^) 

1=3 



1 + av ii>2 — of i 
1 — at> it>2 
1 







(4) 
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is realizable for any a,v\,v 2 G k[x\, . . . ,x m ]. Let "— indicate that we are 
applying elementary operations, and consider the following: 



1 + av ± v 2 



A 







av±v 2 
1 




/ 1 + av\V 2 
av 2 

\ o 

1 V! 

1 v 2 
-av 2 avi 1 




Vi 

v 2 



av\ 1 + av\v 2 





(1 





°1 


/l 





0\ 


)- 





1 


v 2 -»• 





1 







Vo 





1 / 


Vo 





1/ 



(5) 



y af i 1 + auiv 2 
Keeping track of all the elementary operations involved, we get 
A = E 13 (-v 1 )E 23 (-v 2 )E 31 (-av 2 )E 32 (av 1 )E 13 (v 1 )E 23 (v 2 )E 31 (av 2 )E 32 (-av 1 ). (6) 
In general (i.e., for arbitrary i < j), 



B 



I + a 



(0,. .. ,0,^,0,. ..,0,-^,0,...,0) 



\Vn) 



(Here, Vj occurs at the i-th position and — Vi occurs at the 
j-th position.) 

/ 1 ■ • • aviVj ■ ■ ■ —aviVi ■ ■ ■ \ 

: : 

1 + aviVj 



-avi 



V n Vj 



1 — aviVj 



-V n Vi 



1/ 
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(\ 







1 + aViV 



V 



av] 







-av; 



1 — aViVj 







J| Ei^aviV^Eiji—aviVi) 

= E it (-Vi) E jt ( -f j ) (-auj ) (a^ ) E it (vi)E jt (vj ) £ ti (aVj ) ^ (-auj) 
JJ Eufaviv^Ei^-aviVi). (7) 

In the above, t G {1, . . . ,n} can be chosen to be any number other than i 
and j. □ 

Since a Cohn-type matrix is realizable, any product of Cohn-type matrices 
is also realizable. This observation motivates the following generalization of 
the above lemma. 

Definition 3 Let R be a ring and v = (vi, . . . , v n Y G R n for some n G IV. 
Then v is called a unimodular column vector if its components generate R, 
i.e. if there exist gi, . . . , g n G R such that v±gi + • ■ • + v n g n = 1. 



Corollary 1 Suppose that A G SL n (k[xi, . . . , x m ]) with n > 3 can be written 
in the form A = I + v-w for a unimodular column vector v and a row vector 
w over k[xi, . . . , x m ] such that w • v = 0. Then A is realizable. 

Proof: Since v = (t> 1? . . . , v n Y is unimodular, we can find gi,...,g n G 
k[xi, . . . ,x m ] such that v^gi + • • • + v n g n = 1. We can use the effective 
Nullstellensatz to explicitly find these g^s (See ||). This combined with 
w • v = wiVi + • • ■ + w n v n = yields a new expression for w: 

w = J2 a ij( v i e i -^ e j) (8) 

i<j 
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where a^- = Wigj — Wjgi. Now, 

a = Hi 1 + v ■ a n( v j e i- v i e i)) ■ ( 9 ) 

i<j 

Each component on the right hand side of this equation is a Cohn-type matrix 
and thus realizable, so A is also realizable. □ 

Corollary 2 BE i j(a)B~ 1 is realizable for any B G GL n (fc[xi,...,i m ]) 
with n > 3 and a E k[x\, . . . , x m ] . 

Proof: Note that i ^ j, and 

BEij(a)B L = I + (i-th column vector of 5) • a • (j-th row vector of B^ 1 ). 

Let v be the i-th column vector of B and w be a times the j-th row vector 
of B" 1 . Then (z'-th row vector of -B -1 ) • v = 1 implies v is unimodular, and 
w ■ v is clearly zero since % ^ j. Therefore, BE i j(a)B^ 1 — I + v • w satisfies 
the condition of the above corollary, and is thus realizable. □ 



Remark 2 One important consequence of this corollary is that for n > 
3, E n (k[xi, . . . ,x m ]) is a normal subgroup of SL n (k[xi, . . . , x m }), i.e. if 
A G SL n (k[xi, . . . , x m ]) and E e E n (k[xi, . . . , x m ]), then the above corol- 
lary gives us an algorithm for finding elementary matrices E±, . . . , E t such 
that A~ X EA = E X ---E t . 

3 Glueing of Local Realizability 

Let R = k[x±, . . . ,x m -i], X = x m and M G Max(i?) ={ maximal ideals of 
R}. For A G SL n (R[X}), we let A M G SL n (R M [X]) be its image under the 
canonical mapping SL n (R[X]) —> SL h (Rm[X]). Also, by induction, we may 
assume SL n (R) = E n (R) for n > 3. Now consider the following analogue of 
Quillen's theorem for elementary matrices; 

Suppose n > 3 and A G SL n (R[X]). Then A is realizable over 
R[X] if and only if A M G SL n (R M [X]) is realizable over R M [X] 
for every M G Max(i?) . 
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While a non-constructive proof of this assertion is given in |I2| and a more 
general functorial treatment of this Quillen Induction Process can be found 
in |J, we will attempt to give a constructive proof for it here. Since the 
necessity of the condition is clear, we have to prove the following; 

Theorem 1 (Quillen Induction Algorithm) For any given A G SL n (R[X]), 
if A M G E n (R M [X]) for every M G Max(R), then A G E n (R[X]). 



Remark 3 In view of this theorem, for any given A G SL n (R[X}) , now it's 
enough to have a realization algorithm for each Am over Rm[X]. 

Proof: Let a x = (0, . . . , 0) G /c m ~\ and M 1 = {g G k[x x , . . . , x m _i] | #(ai) = 
0} be the corresponding maximal ideal. Then by the condition of the theo- 
rem, Am x is realizable over Rm^X]. Hence, we can write 

^.-n^(|) (io) 

where Cj, dj G R, dj G" M\. Letting r% = ]Jj dj ^ M 1; we can rewrite this as 
A Ml = ( ° 3 ' n * i #3 ' dfc ) e MRr.) C E n (R Ml ). (11) 

Denote an algebraic closure of k by k. Inductively, let aj G k m ~ l be a 
common zero of r 1; . . . , r 3 -_i and Mj = {g G . . . , x m _i] | ^(aj) = 0} be 
the corresponding maximal ideal of R for each j > 2. Define r,- ^ Mj in the 
same way as in the above so that 

A M .eE n (R r .[X]). (12) 

Since aj is a common zero of r±, . . . , r 3 -_i in this construction, we immediately 
see ri, . . . , rj-i G Mj = {g E R \ g(aj) = 0}. But noting rj ^ M,-, we 
conclude that r 3 - ^ rii? + ■ • • + Tj-iR. Now, since the Noetherian condition 
on i? guarantees that we will get to some L after a finite number of steps such 
that riR + • • • + r^R = R, we can use the usual Ideal Membership Algorithm 
to determine when 1# is in the ideal r%R + • • • + r^R. 

Let / be a large natural number (It will soon be clear what large means). 
Then since r[R + • • ■ + r l L R = R, we can use the effective Nullstellensatz 
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to find <7i, . . . , <7£ G -R such that r\g\ + • • • + r l L gL = 1- Now, we express 
A(X) G SX„(.R[X]) in the following way: 

A(X) = A{X-XT\ gi ).[A-\X-Xr[ gi )A(X)] 

= A(X - Xr[ gi - Xr{g 2 ) ■ [A~\X - Xr\ 9l - Xr l 2 g 2 )A(X - Xr[ 9l )] 
\A-\X - Xr[ gi )A{X)\ 

= A(X -jrXrigA • [A^X - £Xr^)A(X - J^Xr^)] • • • 

i=l i=l i=l 

■ ■■[A-\X-Xr[ gi )A{X)]. (13) 

Note here that the first matrix A(X — Y^i=\ ^ r i9i) = A(0) on the right 
hand side is in SL n (R) = E n {R) by the induction hypothesis. What will be 
shown now is that for a sufficiently large /, each expression in the brackets 
in the above equation for A is actually in E n (R[X\), so that A itself is in 
E n (R[X]). To this end, by letting A M . = A { and identifying A G SL n (R[X]) 
with Ai G SL n (R M .[X]), note that each expression in the brackets is in the 
following form: 

Ar\cX)M(c + r l t g)X). (14) 



*Claim: For any c, g G R, we can find a sufficiently large I such that 

A^ 1 (cX)A i ((c + r\g)X) G E n {R[X\) for alH = 1, . . . , L. 

Let 

Di(X, Y, Z) = A'\Y ■ X)M(Y + Z) • X) G E n {R n [X, Y, Z}) (15) 
and write D-i in the form 

D l = \{E Sjtj {b ] + Zf J ) (16) 

3=1 

where bj G R n [X, Y] and fj G R n [X, Y, Z] . From now on, the elementary 
matrix E Sjtj (a) will be simply denoted as E^{a) for notational convenience. 
Now define C p by 

C p =f{E3{b,)eE n {R n [X,Y]). (17) 
j'=i 
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Then the C p 's satisfy the following recursive relations; 

E\h) = d 

E?(b p ) = C-1 X C V (2<p<h) 
C h = I. 

Hence, using E^a + b) = E i j(a)E i j(b), 



(18) 



3=1 

= f[ E j (bj)E j (Zfj) 

3=1 

= [E^hWZffll&ibd&iZh)] ••• [E h (b h )E h (Zf h 
= [C^iZf^C^C^iZh)] ■■■ [C h \C h E h (Zf h )] 

= Uc^izf^c- 1 . 

3=1 



(19) 



Now in the same way as in the proof of Corollary 1 and Corollary 2 of 

section 2, we can write CjE 3 \Zf^)C~ x as a product of Cohn-type matrices, 



i.e. for any given j G {1, . . . , h}, let v = 



of Cj. Then 



\v n ) 



be the Sj-ih column vector 



CjE,. tj (Zfj)Cj 1 = 1] [I + -vZf r < h s{v y e s -v s e y )] (20) 

!<7<<5<n 



for some G R n [X, Y] . Also we can find a natural number I such that 



7 r\ ' 



_ f — 11 



(21) 



for some i/ , a' y5 G R[X, Y], /j G -R[X, Y, Z\. Now, replacing Z by rf g, we see 
that all the Cohn-type matrices in the above expression for CjE\Zfj)Cj l 
have denominator-free entries. Therefore, 



CjEttfgMC^eEniRfrY]). 



(22) 
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Since this is true for each j, we conclude that for a sufficiently large /, 



D t (X,Y,rlg) = JJ C 3 E^r\gf 3 )Cj l G E n (R[X,Y]). 

3=1 

Now, letting Y = c proves the claim. 



(23) 



□ 



4 Reduction to SL^(k[xi, . . . , x m 

Let A G SL n (k[xi, . . . ,x m ]) with n > 3, and v be its last column vec- 
tor. Then v is unimodular. (Recall that the cofactor expansion along 
the last column gives a required relation.) Now, if we can reduce v to 
e n = (0,0,..., 0, 1)* by applying elementary operations, i.e. if we can find 
B G E n (k[x 1: . . . , x m ]) such that B\ = e n , then 



BA 



A 



0\ 







(24) 



\Pl ... Vn-\ 1/ 



for some A e S , L n _i(/c[a;i, . . . , x m }) and p« G . . . , x m ] for i = 1, . . . , n — \. 
Hence, 



BAE nX (--p\) ■ ■ --Bn(n-l)(-Pn-l) 



A 
1 



(25) 



Therefore our problem of expressing A G SX n (£;[:ri, . . . , x m ]) as a prod- 
uct of elementary matrices is now reduced to the same problem for A G 
SL n _i(k[xi, . . . , x m ]). By repeating this process, we get to the problem of 
( V 1 \ 

expressing A = r s G S , L 3 (A;[xi, . . . , x m ]) as a product of elementary 

\o o i) 

matrices, which is the subject of the next section. In this section, we will 
develop an algorithm for finding elementary operations that reduce a given 
unimodular column vector v G (k[xi, . . . ,x m ]) n to e n . Also, as a corollary 
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to this Elementary Column Property, we give an algorithmic proof of the 
Unimodular Column Property which states that for any given unimodular 
column vector v G {k[x\, . . . , x m ]) n , there exists a unimodular matrix B, i.e. 
a matrix of constant determinant, over k[x\, . . . , x m ] such that Bv = e n . 
Lately, A. Logar, B. Sturmfels in [|[ and N. Fitchas, A. Galligo in 0, @ 
have given different algorithmic proofs of this Unimodular Column Property, 
thereby giving algorithmic proofs of the Quillen-Suslin theorem. Therefore, 
our algorithm gives another constructive proof of the Quillen-Suslin theorem. 
The second author has given a different algorithmic proof of the Elementary 
Column Property based on a localization and patching process in [|14[ . 

Definition 4 For a ring R, XJm n (R) = {n- dimensional unimodular column 
vectors over R}. 

Remark 4 Note that the groups GL n (k[x\, . . . , x m }) and E n (k[xi, . . . , x m ]) 
act on the set XJm n (k[xi, . . . ,x m }) by matrix multiplication. 



Theorem 2 (Elementary Column Property) Forn > 3, the group E n (k[xi, . . . , 
acts transitively on the set Um n (k[xi, . . . ,x m }). 



Remark 5 According to this theorem, if\, v' are n-dimensional unimodular 
column vectors over k[xi, . . . ,x m ], then we can find B G E n {k[xi, . . . ,x m }) 
such that Bv = V . Letting v' = e n gives a desired algorithm. 



Corollary 3 (Unimodular Column Property) Forn > 2, the group GL n (k[xi, 
acts transitively on the set Um n (k[xi, . . . , x m ]). 

Proof: For n > 3, the Elementary Column Property cleary implies the Uni- 
modular Column Property since a product of elementary matrices is always 
unimodular, i.e. has a constant determinant. 

If n = 2, for any v = (vi,v 2 Y G Um 2 (A;[xi, . . . , x m ]), find gi,g 2 G 
k[x\, . . . , x m ] such that v \g\ + v 2 g 2 = 1- Then the unimodular matrix U v = 

f V<2 Vl | satisfies U v ■ v = e 2 . Therefore we see that, for any v, w G 
V 9i 92 J 

Um 2 (fc[ii, . . . , x m }), U~ l U v • v = w where U^U V G GL 2 (k[xi, . . . , x m \). □ 
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Let R = k[x\, . . . ,x m -i] and X = x m . Then k[xi,...,x m ] = R[X]. By 

identifying A G SL 2 (R[X\) with ( „ r ^ ) G SL n (R[X]), we can regard 

SX 2 (-R[X]) as a subgroup of SL n (R[X]). Now consider the following theo- 
rem. 



Theorem 3 Suppose v(X) 



G Um n (i?[X]) ; andfi(X) is monic 



\Vn(X)J 

in X. Then there exists B\ G SL 2 (R[X]) and B 2 G E n (R[X]) such that 

//;/;, -V(.V) V(()!. 

Proof: Later □ 

We will use this theorem to prove the Theorem [2L now. 

Proof of Theorem |2|: Since the Euclidean division algorithm for k[x\] 
proves the theorem for m = 1 case, by induction, we may assume the 
statement of the theorem for R = k[xi, . . . ,x m _i]. Let X = x m and v = 

( Vl \ 

■ G Um n (i?[X]). We may also assume that V\ is monic by applying a 
\v n ) 

change of variables (as in the well-known proof of the Noether Normalization 
Lemma). Now by the above Theorem ^|, we can find B~i G SL>2(R[X\) and 
B 2 G E n {R[X\) such that 

B X B 2 ■ v(X) = v(0) G R. (26) 

And then by the inductive hypothesis, we can find B' G E n (R) such that 

B' ■ v(0) = e„. (27) 

Therefore, we get 

v = B^B^B'^en. (28) 
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By the normality of E n (R[X]) in SL n (R[X]) (Corollary 2), we can write 
B^B'- 1 = B"B- V X for some B" G E n (R[X}). Since 



-1 



( p q 
r s 




0\ 




<n-2 







(29) 



for some p,q,r,s G -R[X], we have 

v = B^B^B'- 1 ^ 
= (B 2 l B")B^e n 



(B^B») 



( p 


q ... 


0\ 


(°\ 


r 


s ... 




















: In-2 







U 





/ 


w 



(B^B") 



(30) 



where B 2 B" G £J n (i2[X]). Since we have this relationship for any v G 
Um n (_R[X]), we get the desired transitivity. □ 

Now, we need one lemma to construct an algorithm for the Theorem [3|. 

Lemma 3 Let f\, f2,b,d G R[X] and r be the resultant of fx and f%. Then 
there exists B G SL 2 (R[X]) such that 



B 



A(b) 
Mb) 



fi(b + rd) 
f2(b + rd) 



(31) 



Proof: By the property of the resultant of two polynomials, we can find 
9i, 92 G R[X] such that figi + f 2 g2 = r. Also let s\, s 2 , fa, t 2 G R[X, Y, Z] be 
the polynomials defined by 



ft(X + YZ) 
f 2 (X + YZ) 
9l (X + YZ) 
g 2 (X + YZ) 



f 1 (X)+Ys l (X,Y,Z) 
f 2 (X)+Ys 2 (X,Y,Z) 
9l {X)+Yfa(X,Y,Z) 
g2(X) + Yt 2 (X,Y,Z). 



(32) 
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Now, let 



B n = l+ Sl (b,r,d)-gi(b)+h{b,r,d)- f 2 (b) 

Bn = sx{b,r,d) • g 2 {b) -t 2 (b,r,d) - f x {b) 

B 2 i = s 2 {b,r,d) • g x (b) -ti(b,r,d) ■ f 2 {b) 

B22 = l + s 2 (b,r,d)-g 2 (b)+t 1 (b,r,d)- ft(b). 



(33) 



Then one checks easily that B 
and that B G SL 2 (R[X}). 



Bn B12 
B 2 \ B22 



satisfies the desired property 

□ 



Proof of Theorem g: Let a x = (0, . ..,0) G k m ~ x . Define M x = {g G 
k[xi, . . . ,x m -i] I ^(ai) = 0} and k% = R/M\ as the corresponding maximal 
ideal and residue field, respectively. Since v e (_R[X]) n is a unimodular 
column vector, its image v in (fcifX])" - = ((R/Mi)[X]) n is also unimodular. 
Since k\[X\ is a principal ideal ring, the minimal Grobner basis of its ideal 
< v 2 ,...,v n > consists of a single element, G\. Then v\ and G\ generate 
the unit ideal in k\[X] since V\, v 2 , . . . , v n generate the unit ideal. Using the 
Euclidean division algorithm for /ci[X], we can find E\ G E n -x(ki\X\) such 
that 



\v n J 



( G x \ 


V J 



(34) 



By identifying k\ with a subring of R, we may regard E\ to be an element of 
E n (R[X]) and G\ to be an element of i?[X]. Then, 



1 

Ex 



t Vx \ 

Gx + qx2 

913 



(35) 



for some q 12 , 



V Qln J 

, q ln G Mx[X]. Now, define r x G R by 

rx = Res(vx, Gx + qxz) 

= the resultant of V\ and Gx + qx2 



(36) 
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and find fi,h± G R[X] such that 

fi-v 1 + h 1 -(G 1 + q 12 )=r 1 . (37) 

Since v\ is monic, and v\ and G\ G k\[X] generate the unit ideal, we have 

fi = Res(fi,Gi + gi 2 ) 
= Res^i, Gi) 

^ 0. (38) 

Therefore, r x ^ M x . Denote an algebraic closure of k by k. Inductively, let 
aj G k" 1 ' 1 be a common zero of r±, . . . , 7j_i and M,- be the corresponding 
maximal ideal of i? for each j > 2. Define 7j ^ Mj in the same way as in 
the above. Define also, Ej G E n _i(kj[X]), Gj G kj[X], fj, hj G -R[X], and 
qj 2 , . ■ . , G Mj[X] in an analogous way. Since we let aj be a common zero 
of ri, ...,rj_i in this construction, we see ri,...,rj_i G Mj = {g G i? | 
(?(aj) = 0}. But noting r,- ^ M,-, we conclude that r,- ^ ri_R + • • • + Tj-iR. 
Now, since i? is Noetherian, after a finite number of steps, we will get to some 
L such that r±R + ■ ■ ■ + r^R = R. We can use the effective Nullstellensatz to 
explicitly find those g^s in R such that T\g\ + • • • + r^gi = 1. Define, now, 
b , bi, . . . , bi, G R[X] in the following way: 

b = 

bi = r\g x X 

b 2 = r x g x X + r 2 g 2 X 

b L = r l9l X + r 2 g 2 X + ■ ■ • + r L g L X = X. (39) 
Then these bi 's satisfy the recursive relations: 
b = 

h = h-i+TigiX for i = 1,...,L. (40) 

*Claim: For each i G there exists 5i G SX 2 (-R[^]) and 5,' G 

f? n (i2[A"]) such that v(6i) = S^'v^-i)- 
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If this claim is true, then using E n (R[X\) ■ SL 2 (R[X\) C SL 2 {R[X\) 
E n {R[X\) (Normality of E n (R[X]); Corollary 2), we inductively get 

v(X) = v(b L ) 

= B L B' L v(b L ^) 



= BB>v(b ) 

= BB'v(O) (41) 

for some B G SL 2 {R[X]) and B' G £ n (i?[X]). Therefore it's enough to prove 
the above claim. For this purpose, let Gi = Gi + q^.. Then 

/ <X) \ 



1 

E(X) 



v(X) 



Gi(X) 

\qin(x) j 



(42) 



For 3 < / < n, we have 

qu{h) - qu{h-i) G (bi - 6i_i) • i?[X] 

= r^X-tffX]. (43) 

Since G i? doesn't depend on X, we have 

n = f i (x) Vl (x) + h i (x)G i (x) 

= fi{h-i)v\{bi-i) + ^(6 i _i)G i (6 i _i) 

= a linear combination of and Gj(6j_i) over i?[X]. (44) 

Therefore, we see that for 3 < I < n, 

qaipi) = quipi^i) + a linear combination of t>i(6j_i) and Gj(6j_i) overi?[X]. 
Hence we can find C G E n (R[X\) such that 

( Mh-i) \ 

' 1 o 



Gi(bi-i) 
Qisipi-i] 



\ qin{bi-l) J 
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( Mh-i) \ 
G<(&i-i) 



V qin{k) ) 

Now, by the Lemma |3|, we can find B G SL,2(R[X]) such that 

\Giibi-x) J {Gi(bi)J- 
Finally, define B G SL n (R[X]) as follows: 

Then this B satisfies 

Sv(6i_i) = v(6 i ), 
and by using the normality of J5 n (i2[X]) again, we see that 

b g 5L 2 (it;[x])s n (it:[x]) 

and this proves the claim. 



(45) 



(46) 



(47) 

(48) 

(49) 
□ 



5 Realization Algorithm for SL^(R[X\) 

Now, we want to find a realization algorithm for the matrices of the spe- 

/ p q 

cial type in SL 3 (k[xi, . . . , x m )), i.e. matrices of the form r s 

\0 1 

SL 3 (k[xi, . . . , x m ]). Again, by applying a change of variables, we may as- 
sume that p G k[xi, . . . ,x m ] is a monic polynomial in the last variable x m . 
In view of the Quillen Induction Algorithm developed in the section 3, we 
see that it's enough to develop a realization algorithm for the matrices of 
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/ p q 0\ 

the form r s G SX 3 (i£[Jf]), where i? is now a commutative local 
\0 1) 

ring and p G R[X] is a monic polynomial. A realization algorithm for this 
case was obtained by M.P. Murthy, and we present in the below a slightly 
modified version of the Lemma 3.6 in J3J Suslin's Work on Linear Groups 
over Polynomial Rings and Serre Problem by S.K. Gupta and M.P. Murthy. 



Lemma 4 Let L be a commutative ring, and a,a',b G L. Then, the follow- 
ing s are true. 

1. (a,b) and (a',b) are unimodular over L if and only if (aa',b) is uni- 
modular over L. 

2. For any c,d G L such that aa'd — be = 1, there exist c\, C2, d\, d 2 G L 
such that ad\ — bc\ — 1, a'd 2 — bc^ — 1, and 



I aa' b 

c d | = 
V 1 



fa b 0\ fa' b 

ci d x ■ c 2 d 2 | (mod E 3 (L)). 
\ 1 / V 1 



Proof: (1) If (aa',b) is unimodular over L, there exist hi,h 2 G L such that 
h\ ■ (aa') + h,2 ■ b = 1. Now (h\d) ■ a + h 2 ■ b = 1 implies (a, 6) is unimodular, 
and (h\a) ■ a' + h 2 ■ b = 1 implies (a', 6) is unimodular. 

Suppose, now, that (a, b) and (a', 6) are unimodular over L. Then, we 
can find hi,h 2 , h[, h' 2 G L such that hia + h 2 b = 1, h[a' + /i' 2 6 = 1. Now, let 
gi = hih' x , g 2 = h' 2 + a'h 2 h[, and consider 

giaa' + = h\h\aa! + (/if, + a!h 2 h!^)b 
= h^a'lhia + /^fr) + h' 2 b 
= h^a' + /i' 2 6 

= 1. (50) 
So we have a desired unimodular relation. 

(2) If c, c? G L satisfy aa'd — be = 1, then (aa',b) is unimodular, which in 
turn implies that (a, 6) and (a', 6) are unimodular. Therefore, we can find 
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ci, di, di, d 2 G 1/ such that adi — bc\ — 1 and a'<i 2 — bc 2 — 1. For example, we 
can let 



Ci = c 2 = c, di = a'd, d 2 



ad. 



(51) 



Now, consider 

(ad b 0\ 

c d 
V 



(ad b 0\ 

E 2 i(cdid 2 — d(c 2 + a'cid 2 )) c 2 + dc 1 d 2 d\d 2 

v o o i; 

E 21 (cd 1 d 2 - d(c 2 + a' Cl d 2 ))E 23 {d 2 - 1)E 32 {1)E 23 (-1) 



a b 



a' 6 



ci di £ 23 (1)£ 32 (-1)£ 23 (1) c 2 d 2 



1 



1 



£ 23 (-l)£ 32 (l)£ 23 (a - l)E 31 (-a'c 1 )E 32 (-d 1 ). (52) 
This explicit expression tells us that 



/ ad 


b 




/ a 


b 


°^ 


fa' 


6 




c 


d 


!)■ 


= Ci 




: ■ 


c 2 


d 2 


!) 


V o 






















□ 



Theorem 4 Suppose (R, M) is a commutative local ring, and A 
SL 3 (R[X]) where p is monic. Then A is realizable over R[X]. 



( p q 

r s 

Vo 1 



Proof: By induction on deg(p). If deg(p) = 0, then p = or 1, and A 
is clearly realizable. Now, suppose deg(p) = d > and deg(g) = Z. Since 
p E R[X] is monic, we can find f,g E R[X] such that 

Q = fp + 9, deg(g) < d. (54) 

Then, 

P Q ~ fp 

AE 12 (-f) = | r s-/r | = [ r s - fr |. (55) 

1 
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Hence we may assume deg(g) < d. Now, we note that either p(0) or q(0) 
is a unit in R, otherwise, we would have p(0)s(0) — g(0)r(0) G M that 
contradicts to ps — qr = p(0)s(0) — g(0)r(0) = 1. Let's consider these two 
cases, separately. 

Case 1: When g(0) is a unit. 

Using the invertibility of <?(0), we have 



AE 21 (-q(0)- 1 p(0)) 



p — q(0) 1 p(0)q q 
r - g(0) _1 p(0)s s 
1 



(56) 



So, we may assume p(0) = 0. Now, write p = Xp' . Then, by the above 
Lemma |], we can find c\,di,C2,d2 G R[X] such that Xd\ — qc\ = 1, p'd2 — 
qc2 = 1 and 



(mod E 3 {R[X})) (57) 



Since deg(p') < d, the second matrix on the right hand side is realizable by 
the induction hypothesis. As for the first one, we may assume that q is a unit 
of R since we can assume deg(g) < deg(X) = 1 and q(0) is a unit. And then 

/ X q 





invertibility of q leads easily to an explicit factorization of 



Cl 

V o 



d\ 






1 



into 



elementary matrices. 
Case 2: When q(0) is not a unit. 

First we claim the following; there exist p', q' G R[X] such that deg(p') < 
l,deg(q') < d and p'p — g'g = 1. To prove this claim, we let r G -R be the 
resultant of p and q. Then, there exist /, g G -R[A] with deg(/) < I, deg(g) < 
d such that fp + gq = r. Since p is monic and p, q G -R[A] generate the unit 
ideal, we see that r ^ M, i.e. r £ A*. Now, letting p' — f/r, q' — —g/r 
shows the claim. Also note that the two relations, p'(0)p(0) — q'(0)q(0) = 1 
and g(0) G M, imply p'(0) ^ M. This means q(0) + p'(0) i s a unit. Now, 
consider the following. 



I p q 
r s 
\0 



0\ 


1/ 



p g 

E 21 {rp'-sq') \ q' p' 
1 



21 



E 21 (rp> - sq')E 12 (-V 



( p + q' q + p' 
q' p' 
V 1 



(58) 



Noting that the last matrix on the right hand side is realizable by the Case 1 

( p q 

since g(0) + p'(0) is a unit and deg(p + q') = d, we see that 



also realizable. 



r s 
V0 1 



is 



□ 



6 Eliminating Redundancies 

When applied to a specific matrix, the algorithm presented in this paper 
will produce a factorization into elementary matrices, but this factorization 
may contain more factors than is necessary. The Steinberg relations from 
algebraic i^-theory provide a method for improving a given factorization by 
eliminating some of the unnecessary factors. The Steinberg relations that 
elementary matrices satisfy are 

1. £y(0) = I 

2. E ij (a)E ij (b)=E ij (a + b) 

3. For % ± I, [Eij^Ej^b)] = E i j(a)Eji(b)E i j(—a)Eji(—b) = E u (ab) 

4. For j ^ I, [E l3 (a),E H (b)} = E ij (a)^ i (6)^(-a)^ i (-6) = E^-ab) 

5. For i^pj^l, [E i:j (a),E lp (b)} = ^(a)^ p (6)^(-a)^ p (-6) = /. 

The first author is in the process of implementing the realization algorithm of 
this paper, together with a Redundancy Elimination Algorithm based on the 
above set of relations, using existing computer algebra systems. As suggested 
in [131, an algorithm of this kind has application in Signal Processing since it 



gives a way of expressing a given multidimensional filter bank as a cascade 
of simpler filter banks. 
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